--------------------------------------------------------------------------------
      name:  <unnamed>
       log:  E:\Seagate Sync\VOL\Personal folder\Documents\Research\Past project
> s\cbld\Posted\121020 Website\Data analysis\cbld-analysis01-app.log
  log type:  text
 opened on:  21 Oct 2012, 15:41:37

. 
. 
. * ****************************************************************************
> ********
. * Appendix: Additional results for 'An Empirical evaluation of aggregate model
>  output'
. * ****************************************************************************
> ********
. 
. * Programme:    cbld-analysis01-app.do
. * Project:              Council coalition building
. * Author:               Frank Haege, Department of Politics and Administration
> , University of Limerick
. * Contact:              frank.haege@ul.ie
. 
. * Description
. *************
. * This do-file generates a table and graphs presenting the results of sensitiv
> ity analyses 
. * of the model's fit to the observed consensus rates in the Council. The senst
> itivity analyses    
. * are based on models that incorporate real-world voting weights and threshold
> s, and in which the 
. * sequence of moves of member states is ordered according to voting weights (f
> rom smallest to largest). 
. 
. 
. * Set up Stata
. version 11

. clear all

. macro drop _all

. set linesize 80

. set more off

. 
. * Generate empty matrix to save model fit statistics
. matrix results = J(18,4,.)

. matrix colnames results = "12" "15" "25" "Total"

. matrix rownames results = Observed Null "Model 1" "PRE 1" "Model 2" "PRE 2" "M
> odel 3" "PRE 3" /*
>         */ "Model 4" "PRE 4" "Model 5" "PRE 5" "Model 6" "PRE 6" "Model 7" "PR
> E 7" "Model 8" "PRE 8"

. matrix list results

results[18,4]
             12     15     25  Total
Observed      .      .      .      .
    Null      .      .      .      .
 Model 1      .      .      .      .
   PRE 1      .      .      .      .
 Model 2      .      .      .      .
   PRE 2      .      .      .      .
 Model 3      .      .      .      .
   PRE 3      .      .      .      .
 Model 4      .      .      .      .
   PRE 4      .      .      .      .
 Model 5      .      .      .      .
   PRE 5      .      .      .      .
 Model 6      .      .      .      .
   PRE 6      .      .      .      .
 Model 7      .      .      .      .
   PRE 7      .      .      .      .
 Model 8      .      .      .      .
   PRE 8      .      .      .      .

. 
. 
. * Fit statistics of models with unweighted votes
. ************************************************
. 
. * Set starting value for matrix row number
. local row = 2

. 
. * Loop through model results based on random (a) and ordered sequence of moves
>  (b)
. foreach i in "a" "b" { 
  2. 
.         * Load oberved consensus rates data
.         use "Experiments\Other data\cbld-management05a.dta", clear
  3.         
.         * Rescale consensus rate variable
.         generate oconsensus = consensus/100
  4.         
.         * Drop redundant variables
.         drop consensus
  5. 
.         * Merge observed consensus rates data with predicted consensus rates d
> ata
.         merge m:m rulemembersize using "Experiments\Experiment02\cbld-manageme
> nt03-exp02`i'.dta"
  6.         drop _merge
  7.         
.         * Drop simulation results for 27 member states (under Nice and Lisbon)
.         drop if year == .
  8. 
.         * Calculate squared deviations
.         generate mdev = (oconsensus - consensus)^2
  9.         generate ndev = (oconsensus - econsensus)^2
 10. 
.         * Generate sum of square
.         egen mss = sum(mdev), by(nostates)
 11.         egen tss = sum(ndev), by(nostates)
 12.         egen tmss = sum(mdev)
 13.         egen ttss = sum(ndev)
 14.         drop oconsensus mdev ndev
 15.         
.         * Generate R-squared
.         generate r2 = (tss - mss) / tss
 16.         generate tr2 = (ttss - tmss) / ttss
 17. 
.         * Collapse data set by number of member states
.                 collapse (mean) r2=r2 tr2=tr2 consensus=consensus se=se lowerc
> i=lowerci upperci=upperci /*
>                 */ econsensus=econsensus ese=ese elowerci=elowerci eupperci=eu
> pperci mss=mss tss=tss /*
>                 */ tmss=tmss ttss=ttss, by(rulemembersize)
 18. 
.         * Label variables
.         label var rulemembersize "Rule-membership period"
 19.         label var r2 "Prop. reduction of error for regime periods"
 20.         label var tr2 "Prop. reduction of error for entire period"
 21.         label var consensus "Expected consensus (model)"
 22.         label var se "Standard error (model)"
 23.         label var lowerci "CI lower bound (model)"
 24.         label var upperci "CI upper bound (model)"
 25.         label var econsensus "Expected consensus (null)"
 26.         label var ese "Standard error (null)"
 27.         label var elowerci "CI lower bound (null)"
 28.         label var eupperci "CI upper bound (null)"      
 29. 
.         * Merge with observed data
.         merge 1:m rulemembersize using "Experiments\Other data\cbld-management
> 05b.dta"
 30.         drop _merge
 31. 
.         * Generate variable for running number
.         generate no = _n
 32.         replace no = no - 3.2 if no > 3
 33. 
.         * Merge with observed data
.         merge m:m rulemembersize using "Experiments\Other data\cbld-management
> 05c.dta"
 34.         drop _merge
 35. 
.         * Update running number variable
.         replace no = _n if no == .
 36.         replace no = no - 6.1 if no > 6
 37. 
.         * Label running number variable
.         label var no "Number of Member States"
 38.         label def nol 1 "12" 2 "15" 3 "25", replace
 39.         label val no nol
 40. 
.         * Rescale variables
.         replace r2 = r2 *100
 41.         replace tr2 = tr2 *100
 42.         generate pconsensus = consensus *100
 43.         generate plowerci = lowerci * 100
 44.         generate pupperci = upperci * 100
 45.         generate peconsensus = econsensus * 100
 46.         generate pelowerci = elowerci * 100
 47.         generate peupperci = eupperci * 100
 48. 
.         * Save R-squared measures in locals
.         sort rulemembersize
 49.         mkmat rulemembersize r2 tr2 mconsensus pconsensus peconsensus, matr
> ix(r2)
 50.         matrix list r2
 51.         local r2_12 = string(round(r2[3,2], 1))
 52.         local r2_15 = string(round(r2[6,2], 1))
 53.         local r2_25 = string(round(r2[9,2], 1)) 
 54.         local tr2 = string(round(r2[3,3], 1))
 55. 
.         * Extract results and save them in matrix
.         matrix results[1,1] = r2[2,4]
 56.         matrix results[1,2] = r2[5,4]
 57.         matrix results[1,3] = r2[8,4]
 58.         matrix results[2,1] = r2[3,6]
 59.         matrix results[2,2] = r2[6,6]
 60.         matrix results[2,3] = r2[9,6]
 61.         local ++row
 62.         matrix results[`row',1] = r2[3,5]
 63.         matrix results[`row',2] = r2[6,5]
 64.         matrix results[`row',3] = r2[9,5]
 65.         local ++row
 66.         matrix results[`row',1] = r2[3,2]
 67.         matrix results[`row',2] = r2[6,2]
 68.         matrix results[`row',3] = r2[9,2]
 69.         matrix results[`row',4] = r2[3,3]
 70.         matrix list results
 71. 
.         * Plot model and null-model estimates against observed data
.         twoway (rcap pupperci plowerci no, /*
>                 */ ylabel(0 20:100) xsize(6) /*
>                 */ ytitle("Consensus Rate (Per Cent)") lwidth(medium) lcolor(g
> rey) msize(medium) /*
>                 */ xscale(range(0 3.5)) xlabel(1 2 3, val)) /*
>                 */ (rcap peupperci pelowerci no, /*
>                 */ lwidth(medium) lcolor(grey) msize(medium)) /*
>                 */ (scatter peconsensus no, msymbol(s) mcolor(white)) /*
>                 */ (scatter pconsensus no, msymbol(d) mcolor(white)) /*
>                 */ (scatter mconsensus no, msymbol(dh) mcolor(white)) /*
>                 */ (scatter consensus1994 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus1995 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus1996 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus1997 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus1998 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus1999 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus2000 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus2001 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus2002 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus2003 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus2004 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus2005 no, msymbol(oh) mcolor(white)) /*
>                 */ (scatter consensus2006 no, msymbol(oh) mcolor(white)) /*
>                 */ , legend(order(- "Predicted:" 4 "Computational Model" 3 "Nu
> ll Model" /*
>                 */ - " " - "Observed:" 5 "Mean" 6 "Individual Values" - " " - 
> "Total {it:PRE}: `tr2'%") /*
>                 */ cols(1) ring(1) position(3) region(lcolor(black))) scale(1.
> 2) /*
>                 */ text(10 .3 "{it:PRE} :") text(10 1 "`r2_12'%") text(10 2 "`
> r2_15'%") text(10 3 "`r2_25'%") /*
>                 */ saving("Data analysis\Graphs\cbld-analysis01-app01`i'.gph",
>  replace)
 72. }
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             2
        from master                         0  (_merge==1)
        from using                          2  (_merge==2)

    matched                                13  (_merge==3)
    -----------------------------------------
(2 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     67.737595     94.801201             .     89.200005
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     95.846931     94.801201             .     84.599998
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     96.698654     94.801201             .     85.199997

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2          .          .          .          .
   PRE 2          .          .          .          .
 Model 3          .          .          .          .
   PRE 3          .          .          .          .
 Model 4          .          .          .          .
   PRE 4          .          .          .          .
 Model 5          .          .          .          .
   PRE 5          .          .          .          .
 Model 6          .          .          .          .
   PRE 6          .          .          .          .
 Model 7          .          .          .          .
   PRE 7          .          .          .          .
 Model 8          .          .          .          .
   PRE 8          .          .          .          .
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app01a.gph saved)
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             2
        from master                         0  (_merge==1)
        from using                          2  (_merge==2)

    matched                                13  (_merge==3)
    -----------------------------------------
(2 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     65.425591     94.499641             .     89.700005
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     95.609192     94.499641             .            84
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     96.577759     94.499641             .     81.599998

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2  89.700005         84  81.599998          .
   PRE 2  65.425591  95.609192  96.577759  94.499641
 Model 3          .          .          .          .
   PRE 3          .          .          .          .
 Model 4          .          .          .          .
   PRE 4          .          .          .          .
 Model 5          .          .          .          .
   PRE 5          .          .          .          .
 Model 6          .          .          .          .
   PRE 6          .          .          .          .
 Model 7          .          .          .          .
   PRE 7          .          .          .          .
 Model 8          .          .          .          .
   PRE 8          .          .          .          .
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app01b.gph saved)

. 
. 
. * Fit statistics of models with weighted votes and several voting thresholds
. ****************************************************************************
. 
. * Loop through model results based on 1 to 3 voting thresholds
. foreach j of num 1 2 3 { 
  2.         * Loop through model results based on random (a) and ordered sequen
> ce of moves (b)
.         foreach i in "a" "b" {
  3. 
.                 * Load observed consensus rates data
.                 use "Experiments\Other data\cbld-management05a.dta", clear
  4.                 
.                 * Rescale consensus rate variable
.                 generate oconsensus = consensus/100
  5.                 
.                 * Drop redundant variables
.                 drop consensus
  6. 
.                 * Merge observed consensus rates data with predicted consensus
>  rates data
.                 merge m:m rulemembersize using "Experiments\Experiment03\cbld-
> management04-exp03`i'`j'.dta"
  7.                 drop _merge
  8.         
.                 * Drop simulation results for 27 member states (under Nice and
>  Lisbon)
.                 drop if year == .
  9.         
.                 * Calculate squared deviations
.                 generate mdev = (oconsensus - consensus)^2
 10.                 generate ndev = (oconsensus - econsensus)^2
 11. 
.                 * Generate sum of square
.                 egen mss = sum(mdev), by(nostates)
 12.                 egen tss = sum(ndev), by(nostates)
 13.                 egen tmss = sum(mdev)
 14.                 egen ttss = sum(ndev)
 15.                 drop oconsensus mdev ndev
 16. 
.                 * Generate R-squared
.                 generate r2 = (tss - mss) / tss
 17.                 generate tr2 = (ttss - tmss) / ttss
 18. 
.                 * Collapse dataset by number of member stats
.                 collapse (mean) r2=r2 tr2=tr2 consensus=consensus se=se lowerc
> i=lowerci upperci=upperci /*
>                         */ econsensus=econsensus ese=ese elowerci=elowerci eup
> perci=eupperci mss=mss tss=tss /*
>                         */ tmss=tmss ttss=ttss, by(rulemembersize)
 19. 
.                 * Label variables
.                 label var rulemembersize "Rule-membership period"
 20.                 label var r2 "Prop. reduction of error for regime periods"
 21.                 label var tr2 "Prop. reduction of error for entire period"
 22.                 label var consensus "Expected consensus (model)"
 23.                 label var se "Standard error (model)"
 24.                 label var lowerci "CI lower bound (model)"
 25.                 label var upperci "CI upper bound (model)"
 26.                 label var econsensus "Expected consensus (null)"
 27.                 label var ese "Standard error (null)"
 28.                 label var elowerci "CI lower bound (null)"
 29.                 label var eupperci "CI upper bound (null)"      
 30. 
.                 * Merge with observed consensus rates data
.                 merge 1:m rulemembersize using "Experiments\Other data\cbld-ma
> nagement05b.dta"
 31.                 drop _merge
 32. 
.                 * Generate variable for running number
.                 generate no = _n
 33.                 replace no = no - 3.2 if no > 3
 34. 
.                 * Merge with observed consensus rates data
.                 merge m:m rulemembersize using "Experiments\Other data\cbld-ma
> nagement05c.dta"
 35.                 drop _merge
 36. 
.                 * Update running number variable
.                 replace no = _n if no == .
 37.                 replace no = no - 6.1 if no >= 7
 38. 
.                 * Label running number variable
.                 label var no "Number of Member States"
 39.                 label def nol 1 "12" 2 "15" 3 "25", replace
 40.                 label val no nol
 41. 
.                 * Rescale variables
.                 replace r2 = r2 *100
 42.                 replace tr2 = tr2 *100
 43.                 generate pconsensus = consensus *100
 44.                 generate plowerci = lowerci * 100
 45.                 generate pupperci = upperci * 100
 46.                 generate peconsensus = econsensus * 100
 47.                 generate pelowerci = elowerci * 100
 48.                 generate peupperci = eupperci * 100
 49. 
.                 * Save R-squared statistics in locals
.                 sort rulemembersize
 50.                 mkmat rulemembersize r2 tr2 mconsensus pconsensus peconsens
> us, matrix(r2)
 51.                 matrix list r2
 52.                 local r2_12 = string(round(r2[3,2], 1))
 53.                 local r2_15 = string(round(r2[6,2], 1))
 54.                 local r2_25 = string(round(r2[9,2], 1)) 
 55.                 local tr2 = string(round(r2[3,3], 1))
 56. 
.                 * Extract results and save them in matrix
.                 matrix results[1,1] = r2[2,4]
 57.                 matrix results[1,2] = r2[5,4]
 58.                 matrix results[1,3] = r2[8,4]
 59.                 matrix results[2,1] = r2[3,6]
 60.                 matrix results[2,2] = r2[6,6]
 61.                 matrix results[2,3] = r2[9,6]
 62.                 local ++row
 63.                 matrix results[`row',1] = r2[3,5]
 64.                 matrix results[`row',2] = r2[6,5]
 65.                 matrix results[`row',3] = r2[9,5]
 66.                 local ++row
 67.                 matrix results[`row',1] = r2[3,2]
 68.                 matrix results[`row',2] = r2[6,2]
 69.                 matrix results[`row',3] = r2[9,2]
 70.                 matrix results[`row',4] = r2[3,3]
 71.                 matrix list results
 72. 
.                 * Plot model and null-model estimates against observed data
.                 twoway (rcap pupperci plowerci no, /*
>                         */ ylabel(0 20:100) xsize(6) /*
>                         */ ytitle("Consensus Rate (Per Cent)") lwidth(medium) 
> lcolor(grey) msize(medium) /*
>                         */ xscale(range(0 3.5)) xlabel(1 2 3, val)) /*
>                         */ (rcap peupperci pelowerci no, /*
>                         */ lwidth(medium) lcolor(grey) msize(medium)) /*
>                         */ (scatter peconsensus no, msymbol(s) mcolor(white)) 
> /*
>                         */ (scatter pconsensus no, msymbol(d) mcolor(white)) /
> *
>                         */ (scatter mconsensus no, msymbol(dh) mcolor(white)) 
> /*
>                         */ (scatter consensus1994 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus1995 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus1996 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus1997 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus1998 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus1999 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus2000 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus2001 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus2002 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus2003 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus2004 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus2005 no, msymbol(oh) mcolor(white
> )) /*
>                         */ (scatter consensus2006 no, msymbol(oh) mcolor(white
> )) /*
>                         */ , legend(order(- "Predicted:" 4 "Computational Mode
> l" 3 "Null Model" /*
>                         */ - " " - "Observed:" 5 "Mean" 6 "Individual Values" 
> - " " - "Total {it:PRE}: `tr2'%") /*
>                         */ cols(1) ring(1) position(3) region(lcolor(black))) 
> scale(1.2) /*
>                         */ text(10 .3 "{it:PRE} :") text(10 1 "`r2_12'%") text
> (10 2 "`r2_15'%") text(10 3 "`r2_25'%") /*
>                         */ saving("Data analysis\Graphs\cbld-analysis01-app02`
> i'`j'.gph", replace)
 73.         }
 74. }
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                13  (_merge==3)
    -----------------------------------------
(0 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     81.337608     94.975685             .     85.799995
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     95.646103     94.975685             .     85.399994
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     95.521072     94.975685             .     88.200005

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2  89.700005         84  81.599998          .
   PRE 2  65.425591  95.609192  96.577759  94.499641
 Model 3  85.799995  85.399994  88.200005          .
   PRE 3  81.337608  95.646103  95.521072  94.975685
 Model 4          .          .          .          .
   PRE 4          .          .          .          .
 Model 5          .          .          .          .
   PRE 5          .          .          .          .
 Model 6          .          .          .          .
   PRE 6          .          .          .          .
 Model 7          .          .          .          .
   PRE 7          .          .          .          .
 Model 8          .          .          .          .
   PRE 8          .          .          .          .
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app02a1.gph saved)
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                13  (_merge==3)
    -----------------------------------------
(0 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     82.694412     96.018875             .     85.399994
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     96.602493     96.018875             .     82.200005
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     96.756142     96.018875             .            86

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2  89.700005         84  81.599998          .
   PRE 2  65.425591  95.609192  96.577759  94.499641
 Model 3  85.799995  85.399994  88.200005          .
   PRE 3  81.337608  95.646103  95.521072  94.975685
 Model 4  85.399994  82.200005         86          .
   PRE 4  82.694412  96.602493  96.756142  96.018875
 Model 5          .          .          .          .
   PRE 5          .          .          .          .
 Model 6          .          .          .          .
   PRE 6          .          .          .          .
 Model 7          .          .          .          .
   PRE 7          .          .          .          .
 Model 8          .          .          .          .
   PRE 8          .          .          .          .
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app02b1.gph saved)
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                13  (_merge==3)
    -----------------------------------------
(0 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     81.337608     94.173126             .     85.799995
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     95.696671     94.173126             .     85.399994
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     92.130264     94.173126             .          92.5

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2  89.700005         84  81.599998          .
   PRE 2  65.425591  95.609192  96.577759  94.499641
 Model 3  85.799995  85.399994  88.200005          .
   PRE 3  81.337608  95.646103  95.521072  94.975685
 Model 4  85.399994  82.200005         86          .
   PRE 4  82.694412  96.602493  96.756142  96.018875
 Model 5  85.799995  85.399994       92.5          .
   PRE 5  81.337608  95.696671  92.130264  94.173126
 Model 6          .          .          .          .
   PRE 6          .          .          .          .
 Model 7          .          .          .          .
   PRE 7          .          .          .          .
 Model 8          .          .          .          .
   PRE 8          .          .          .          .
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app02a2.gph saved)
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                13  (_merge==3)
    -----------------------------------------
(0 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     82.694412     95.220604             .     85.399994
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     96.817711     95.220604             .     82.200005
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     92.911179     95.220604             .     92.099998

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2  89.700005         84  81.599998          .
   PRE 2  65.425591  95.609192  96.577759  94.499641
 Model 3  85.799995  85.399994  88.200005          .
   PRE 3  81.337608  95.646103  95.521072  94.975685
 Model 4  85.399994  82.200005         86          .
   PRE 4  82.694412  96.602493  96.756142  96.018875
 Model 5  85.799995  85.399994       92.5          .
   PRE 5  81.337608  95.696671  92.130264  94.173126
 Model 6  85.399994  82.200005  92.099998          .
   PRE 6  82.694412  96.817711  92.911179  95.220604
 Model 7          .          .          .          .
   PRE 7          .          .          .          .
 Model 8          .          .          .          .
   PRE 8          .          .          .          .
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app02b2.gph saved)
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                13  (_merge==3)
    -----------------------------------------
(0 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     81.337608     93.840393             .     85.799995
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     95.655869     93.840393             .     85.399994
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     90.901337     93.840393             .     93.599998

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2  89.700005         84  81.599998          .
   PRE 2  65.425591  95.609192  96.577759  94.499641
 Model 3  85.799995  85.399994  88.200005          .
   PRE 3  81.337608  95.646103  95.521072  94.975685
 Model 4  85.399994  82.200005         86          .
   PRE 4  82.694412  96.602493  96.756142  96.018875
 Model 5  85.799995  85.399994       92.5          .
   PRE 5  81.337608  95.696671  92.130264  94.173126
 Model 6  85.399994  82.200005  92.099998          .
   PRE 6  82.694412  96.817711  92.911179  95.220604
 Model 7  85.799995  85.399994  93.599998          .
   PRE 7  81.337608  95.655869  90.901337  93.840393
 Model 8          .          .          .          .
   PRE 8          .          .          .          .
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app02a3.gph saved)
nostates was byte now float
rulemembersize was byte now float

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                13  (_merge==3)
    -----------------------------------------
(0 observations deleted)

    Result                           # of obs.
    -----------------------------------------
    not matched                             6
        from master                         3  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)

    Result                           # of obs.
    -----------------------------------------
    not matched                             9
        from master                         6  (_merge==1)
        from using                          3  (_merge==2)

    matched                                 0  (_merge==3)
    -----------------------------------------
(3 real changes made)
(3 real changes made)
(3 real changes made)
(3 real changes made)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)
(6 missing values generated)

r2[9,6]
    rulemember~e            r2           tr2    mconsensus    pconsensus
r1     .80000001             .             .             .             .
r2     .89999998             .             .            75             .
r3             1     82.694412     95.393303             .     85.399994
r4           1.8             .             .             .             .
r5           1.9             .             .     81.666664             .
r6             2     96.827209     95.393303             .     82.200005
r7           2.8             .             .             .             .
r8     2.9000001             .             .            86             .
r9             3     93.582458     95.393303             .     91.399994

     peconsensus
r1             .
r2             .
r3            50
r4             .
r5             .
r6            50
r7             .
r8             .
r9            50

results[18,4]
                 12         15         25      Total
Observed         75  81.666664         86          .
    Null         50         50         50          .
 Model 1  89.200005  84.599998  85.199997          .
   PRE 1  67.737595  95.846931  96.698654  94.801201
 Model 2  89.700005         84  81.599998          .
   PRE 2  65.425591  95.609192  96.577759  94.499641
 Model 3  85.799995  85.399994  88.200005          .
   PRE 3  81.337608  95.646103  95.521072  94.975685
 Model 4  85.399994  82.200005         86          .
   PRE 4  82.694412  96.602493  96.756142  96.018875
 Model 5  85.799995  85.399994       92.5          .
   PRE 5  81.337608  95.696671  92.130264  94.173126
 Model 6  85.399994  82.200005  92.099998          .
   PRE 6  82.694412  96.817711  92.911179  95.220604
 Model 7  85.799995  85.399994  93.599998          .
   PRE 7  81.337608  95.655869  90.901337  93.840393
 Model 8  85.399994  82.200005  91.399994          .
   PRE 8  82.694412  96.827209  93.582458  95.393303
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(note:  named style grey not found in class color, default attributes used)
(file Data analysis\Graphs\cbld-analysis01-app02b3.gph saved)

. 
. * TABLE A1.     Observed vs. predicted consensus rates for various model speci
> fications
. ******************************************************************************
> *****
. matrix list results, format(%2.0f)

results[18,4]
             12     15     25  Total
Observed     75     82     86      .
    Null     50     50     50      .
 Model 1     89     85     85      .
   PRE 1     68     96     97     95
 Model 2     90     84     82      .
   PRE 2     65     96     97     94
 Model 3     86     85     88      .
   PRE 3     81     96     96     95
 Model 4     85     82     86      .
   PRE 4     83     97     97     96
 Model 5     86     85     93      .
   PRE 5     81     96     92     94
 Model 6     85     82     92      .
   PRE 6     83     97     93     95
 Model 7     86     85     94      .
   PRE 7     81     96     91     94
 Model 8     85     82     91      .
   PRE 8     83     97     94     95

. 
. 
. log close
      name:  <unnamed>
       log:  E:\Seagate Sync\VOL\Personal folder\Documents\Research\Past project
> s\cbld\Posted\121020 Website\Data analysis\cbld-analysis01-app.log
  log type:  text
 closed on:  21 Oct 2012, 15:41:51
--------------------------------------------------------------------------------
